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ABSTRACT 



A cylindrical acoustic resonator which has been developed by Dr. 
William J. Toulis of the Signal Propagation Division of the United States 
Navy Electronic Laboratory for measurement of the velocity and absorption 
characteristics of an enclosed fluid medium is discussed. The cavity is 
constructed of compliant tubes, thus presenting a pressure release surface 
at the boundaries. Theory is developed indicating that the resonant fre- 
quency of the cavity is directly proportional to the velocity of sound 
in the medium while the Q of the resonator can be used to determine 
attenuation. 

The writer wishes to express his appreciation for the guidance and 
encouragement of Professor Herman Medwin of the United States Naval Post- 
graduate School. He further wishes to acknowledge the assistance given 
him by Dr. Toulis and other members of the Signal Propagation Division 
of the United States Navy Electronics Laboratory in this investigation. 
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1. Introduction. 

The success of sonic echo»ranglng methods depends, In part, upon the 
velocity and attenuation profiles In the transmitting medium between the 
detecting ship and the target. The development of an Instrument which 
could directly measure these properties of the sea would be of great 
assistance. 

The velocity of sound in sea water under varying physical conditions 
has been investigated both experimentally and theoretically by many per- 
sons. For the experimental determinations, much of the early work was 
pointed toward accurate determination of depth by sonic sounding methods. 
Transmission time between a transmitting and receiving hydrophone was 
measured in an area where salinity and temperature were as uniform as 
possible, and in deep enough water so that bottom reflections were mini- 
mized. In this type of determination, measurement of the distance is 
the limiting factor. 

More recently, a number of velocity meters has been developed for 
this purpose. The most successful for laboratory use has been the inter- 
ferometer, where a standing wave is set up in a cell by utilizing a 
crystal and reflector or two reflectors. This then provides a measure of 
the wave length and hence, the velocity. However, no completely satis- 
factory operational-type system is presently available. ^ The most widely 
used device to determine sound velocity profiles is the bathythermograph, 
which measures temperatures rather than velocities. Its use depends on 
the fact that the velocity of sound is largely determined by the tempera- 
ture of the sea, and only to a lesser extent by salinity and pressure. 
Urick^ obtained good general agreement between velocity-depth plots 
obtained with an interferometer and those computed from bathythermograph 
and thermocouple observations for sea states of zero to three and depths 

Sussman, Study and Evaluation of ONR-Sponsored Sound-Velocity 
Meters, USL Research Report No. 299, 1 March 1956 

^R. J. Drick, An Acoustic Interferometer For the Measurement of 
Sound Velocity In the Ocean, USNR & S Lab. Report No. S-18 
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as great as 150 feet. However, the bathythermograph is Insensitive to 
temperature microstructures, ^ and as such, does not give an accurate 
indication of transmission characteristics.^ 

A great deal of continuing effort has been given to the theoretical 
determination of the velocity of sound from the physical constants of 
the medium as given by the theory of thermodynamics. Much of the 
work is, of necessity, empirical in nature and has been complicated by a 
lack of sufficiently accurate experimental data.^ However, recent im> 
provements in laboratory measuring techniques should rapidly improve this 
situation.^ 

Techniques for measuring attenuation of sound by sea water are even 
less satisfactory. The absorption at low frequencies is sufficiently 
small to make measurement difficult without considering a large expanse 
of sea. Here again, determination of distance is a limiting factor. 
Further isolation of the specific causes of attenuation is next to im> 
possible without a device which can determine the absorption character- 
istics of small samples. 

The velocity meter under discussion depends, for its operation, on 
the fact that the resonant frequency of a water column in a particular 

%rick, ibid. 

^L. Liebermann, The Effect of Temperature Inhomogeneities in the 
Ocean On the Propagation of Sound, J. Ac. Soc. Am., 23, pp 563 (1951) 

^S. Kuwahara, The Velocity of Sound in Sea Water and Calculation 
of the Velocity for Use in Sonic Sounding, Japanese J. Astron. & Geophys., 
Vol. 16, No. 1 (1938) 

^D. J. Matthews, Tables of the Velocity of Sound in Pure Water and 
Sea Water for Use in Echo-Sounding & Sound-Ranging, Hydrographic Dept. 
Admiralty, London, H.D. 282 (1927) 

^V. A. Del Grosso, The Velocity of Sound in Sea Water at Zero 
Depth, NRL Report 4002, II June 1952 

®R. T. Beyer, Formulas for Sound Velocity in Sea Water, J. Mar. 
Research, Vol. 13, pp 113 (1954) 

^V. A. Del Grosso, E. J. Smura and P. T. Fougere, Accuracy of 
Ultrasonic Interferometer Velocity Determinations, NRL Report 4439, 

6 December 1954 



2 



mode of vibration is directly proportional to the velocity of sound in 
the water column. The technique employed is not new in Itself. A hard 
wall cavity resonator was developed by the University of Michigan in 1952. 
However, this device was extremely sensitive to wall cleanliness due to 
the relatively high sonic pressures at the boundaries, characteristic of 
rigid walled devices. Thus, although the Instrument provides an accurate 
measure of the velocity when used in the laboratory, it is very difficult 
to use under operational conditions. 

Use of a pressure release walled cylinder suggests Itself; and, in 
fact, thin walled cylinders have been used for measurement of the atten- 
uation of marine sediments . Recent developments in the field of 
compliant structures by Dr. William J. Toulls of the Signal Propagation 
Division of the United States Navy Electronics Laboratory, San Diego, 
California, enabled him to adapt this technique to the present operational 
device, where, by measuring the resonant frequency and the "Q" of the 
cavity, the desired characteristics of a small sample of water can be 
determined to a high degree of accuracy. 

The purpose of this paper is to discuss the theoretical aspects of 
the pressure release walled cylindrical cavity constructed by Toulls, 
including a somewhat systematic study of the loss problem. 



K. Brown, The Development of an Underwater Sound Velocity Meter 
Using a Cylindrical Resonator, ProJ . 2021, final report, 30 April 1953, 
Eng. Research Inst. Univ. of Michigan 

J. Toulls, Theory of a Resoneince Method to Measure the Acoustic 
Properties of Sediments, NBL Report 676, 5 March 1956 

Shumway, A Resonant Chamber Method for Sound Velocity and 
Attenuation Measurement in Sediments, Geophys. Vol XXI, No. 2, pp 305- 
319, April 1956 



2. General Wave Equation. 

Before proceeding to the actual analysis of the velocity meter ^ It 
Is necessary to develop the equations governing the transmission of 
acoustic waves through a fluid medium. It Is assumed that the medium Is 
locally homogeneous, continuous except at the boundaries, and isotropic 
when at rest. It Is further assumed that the condition of the medlttm can 
be fully described by establishing the functional dependence of pressure, 
vector velocity, density and temperature on time and the space coordin- 
ates. A set of equations consisting of 

a) the continuity equation expressing the conservation of mass 

b) a force equation based on conservation of momentum 

c) a heat exchange equation based on conservation of energy 

will be developed In their complete form to present a rigorous point of 
departure for making the necessary approximations for the study of small- 
signal effects. 

To derive the first of the required relationships, consider an 
elemental volume V bounded by a closed surface The amount of fluid 
crossing the surface In time Is given by the surface Integral: 



During this same period, the qixantlty of fluid originally contained In V 
will have diminished by the amount: 



where the negative sign Is selected because density Is a decreasing 
function of time. Since mass must be conserved, net outflow must equal 
the loss of mass and the two expressions can be equated. The surface 
Integral is then converted by use of the divergence theorem giving: 



Since the Integrand Is continuous and the volume arbitrary, this may be 
written as: 



where 



outflow 



• At< j f V • d ^ 



V ■ particle velocity 
~P 3 point fluid density 



loss 3 -A-t f If dV 
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which is one form of the well-known equation of continuity. 
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Considering the density is s function of both space and time, the 
total derlvitive of with respect to *0 can be written in vector 
notation as: 



Dt/ “ ^ ' 
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where 



.D,„. 



+■ V' g 






utilising the notation of Stokes 
rivitlve and is defined by 

p - ^ 

Dt> s t 

Substituting from equation 2.1 into 2.2 gives 

+• -f V V = o 

D b 

which is Euler's fora of the continuity equation. 



2.2 

denotes the material de> 



). 



2.3 



To develop the second of the basic relationships, consider Cauchy's 
14 

first law of motion: 



where : 






^Db/ 






2.4 



px a extraneous force vector, a function of position and 
time 

Va. » velocity component of 
stress tensor. 

For an isotropic mediisn where stress is a linear function of the rate of 
deformation, the stress tensor can be resolved into a sum of a scalar 
and a residual stress tensor as follows: 



- f>&A^ 4- Ua^ (l3A^=t.^A) 2.5 

where 

= pressure 

LAa^m residual stress 
Sa^ b Kronecker delta 

hydrostatic tension. 



^^C. Truesdall, The Mechanical Foundations of Elasticity and Fluid 
Dynamics, J. Rat. Mech. &Anal., 1, pp 125, (1952) 

^^F. V. Hunt, Notes on the Exact Equations Governing The Propaga- 
tion of Sound in Fluids, J. Acous. Soc. of Am., 27, pp 1022, Nov. 1955. 
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To a first approximation, the residual stress can be expressed by 
the linear terms of a power expansion In the viscosity coefficients as: 

X 2.6 

where : 

c(a^ s rate of deformation tensor 
^ s dynamic shear viscosity coefficient 
X B dllatatlonal viscosity coefficient 
and cl 1> defined by 

Combining the equations 2.S, 2.6, and 2.7 and substituting into 2.4 yields 
the vector force equation: 

-f (^* v) V vf + 1X+ a./j)v(v>v) 

~/;VX(vX V) + (X7-'V)VX -V 2-(V/J'V)V 2.8 

+ V/OK (vx \7). 

To develop the third of the fundamental equations. It Is necessary 
to establish a series of relationships based on the law of conservation 
of energy. The total energy, , of a given system can be expressed as 
the sum of the kinetic energy, H. Internal energy, £x« 

Ex=r-P6dv 

'X) - specific volume « '/ -f 
0, • Internal energy density 
Y * volume of the system. 

The material rate of change of the total energy Is then the sum of 
the rates at which work is done upon the system; that Is: 

K Ex) ‘dS 2.10 

where 

= heat flux vector. 

This postulates that the rate of doing mechanical work can be divided 
Into two portions: that arising from the extraneous force and that arls* 

Ing from the forces of material continuity expressed In terms of an 
equivalent stress vector acting on the surface ?. The third term ex- 
presses the rate at which thermal energy enters or leaves the volume. 
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Applying the divergence theorem end utilizing equation 2.4 reduces 2.10 
to 

~ 2.11 

This is the Fourier-Kirchho££-C. Neumann energy equation. Using the 
relationships from equations 2.5 and 2.7, 2.11 can be written as 

-f -f- UA^dA>)r 2.12 

d)(A l! J 

Considering the internal energy density, S , to be a function of S snd 
U gives 




From the first law of thermodynamics. 

Thus, equation 2.13 can be written as 

-p£i. = _ ^V-V 2.14 

Defining , a viscous dissipation function as 



cl 

and eliminating f£| between equations 2.12 and 2.14 gives: , 

= d)^ - ^3 ^ 2.15 

Next, the effects of radiation emd conduction must be considered. 
Since the divergence of gives the energy transferred away from the 

volume element, it must account for energy loss by either conduction or 
radiation. 

For conduction, the Fourier relation is used: 

where 




- thermal conductivity. 

Taking the divergence of each side gives: 

_ /y' v^T “ (v T)* (v/f ) 

b X.A ' 

For the radiation effects, Newton's law of cooling states: 



^ fCv^(T-Ti) 



2.16 



2.17 



7 



where 



^ • radiation coefficient Introduced by Stolces, Phil. Mag. 

(4) 1, pp 305-317 (1851) 

Cv • specific heat at constant volume s ^ V 

(T - Ti) - local temperature excess. 

Before proceeding, the term appearing in equation 2.15 must 

be expressed in terms of , ‘X) , and 'J' . Following the methods of 

Hunt^^, this is given by 









DT 



2.18 



I 

nr Op / cv 
Cp= 

j2> * coefficient of termal expansion 

= -fOV>T)i= 

Combining equation 2.18 with 2.15, 2.16 and 2.17 gives the third of the 
required equations: 



-fCv 



DT 

Dti 



-h -f Cn 



nr— I 
13 



VV-/<V T 



2.19 



-VT-V/r + fCv^ (r- To) - 4>n = O 



^^Hunt, op. cit., pp 1023. 
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3. Small Signal Acoustic Equations. 

The small signal acoustic equations are obtained by arbitrarily 
"linearising" the general eqtiatlons of motion. For this purpose, the 
dependent variables will be replaced by the sum of their equilibrium or 
zero-order values and the first-order variational components. Thus, 

■P= -fo + <^= 

T=T+-©, v=o + ^ 

On the assumption that and constants and that '\/ 

is itself a first-order term (i.e.'^^=r O ) then: 

vf = Vf, = v-ip, 

vT=v9, v-v = '^■v, 

Substituting these relations Into equations 2.1, 2.8 and 2.19 and re- 
taining only the first-order terms gives: 

+ "fo (V-Vi) = o 3.1 



+ V'p, -(yUtvi)v(V-V,) +yUVX(VX.V,) =: Q 



3.2 



3.3 



- -foCv^ O, “ 

where the extraneous force has been omitted since It primarily affects 
the equilibrium configuration and 

N= a ) the viscosity number. 

Considering /p as a function of -p and ""P gives: 



= (If 



since: 



and 



(M!\ 

\ 0~fy rp 



X 



^ _pR 

} t )^ IdfAiT/f 



equation 3.4 becomes 






T 



v-p, i- -fo ve, 
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and 3.2 can be rewritten as 

-^Nv(v-v,) 3_3 

+/J v/(v ^ V,"^ =0. 

These equations (3.1, 3.3, 3.5) take Into account the first-order effects 
of shear and dllatatlonal viscosity and of heat conduction and their 
solution Is discussed In Section 5. They can be reduced to the classical 
small-signal wave equation by neglecting the effects of heat conduction 
and viscosity m /j m)\ mQ) and assuming adiabatic behaviour of 

the medium giving 



^ = 



Ocf- 



3.6 
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4. Solution of the Simple Wave Equation. 

Before considering the more complex equations derived in the pre» 
caeding section, a solution will first be obtained for the small-signal 
wave equation. The results obtained will serve as a valuable guide when 
interpreting the solutions obtaioed in later sections. 

The simple wave equation is given by; 



7. 2. 

- C V -p 



Where V is the Laplacian operator and O is the free space 
velocity. Assuming a time dependence of the form c where 
equation 4.1 takes the form of a scalar Helmholtz equation; 



4.1 

sound 

cr = C0 1- Ao(.^ 



where 



K 



- 

“ o 



_ CO 4- ACX ___ 



4.2 



O 



K 



/\0< 



Using the method of separation of variables with the Laplacian operator 
in cylindrical coordinates gives; 

Jl” 4- -L ^ ± ja'‘ 4- 4. K ^ ~ O A . 

where denotes denotes , etc. By the standard argument 

for the method of separation of variables, equation 4.3 may be rewritten 



as 



R 



± K 

r ~tr 



w 



+ _L -i- -V- = 

® t- t\ 



This then produces two equations ; 

+-k4- + - K: 



'Z. 

'L 



= K. 






T' R 



i" 









4.4a 



4.4b 



where equation 4.4b is the familiar differential equation of simple 
harmonic motion with the solution: 

In general, Ke considered complex; that is 

Ke= ^*^2. Kar.' 

Thus, equation 4.5 may be rewritten as 

^ -f Ke ^ 4.5a 
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Equation 4.4a may be separated utilizing the same techniques giving: 

= O 4.6a 

-|.V O 4.6b 

where solutions to 4.6 may be written in the form: 

^ ~ Cm COS(rr.e+ ^') ^-7 

Since @ must be single valued in 0 , VT\ must be an Integer. 

Equation 4.6b is recognized as a form of Bessel's equation providing 
a solution (with tTi as an integer) of the form: 

■R= JmCKT"^) + El Ny^(Kv'r) 4.8 

where 3m(^vn^® ® Bessel function of the first kind and K1 vk(KyV ® 
Bessel function of the second kind*’®, and iXy — V\ “tY? • The boundary 
conditions which roust be satisfied by equations 4.5a, 4.7 and 4.8 are: 

a. The pressure is finite throughout the interior of the cylinder. 

b. The walls are assumed to be pressure release surfaces requiring 
that J|p=:0 at the boundaries of the cylinder. 

Selecting a coordinate system with Z-O at the bottom of the cylinder, 
and for a cylinder of radius a and length ^ , the following conditions 
apply: 

a. Coefficient of the cosine term in equation 4.5a must be zero 
since ^ is zero at Z.— O* 

b. Periodicity in Z. is given by where a 1,2,3- since /p 

is zero at Z ■= L- . 

e. The coefficient of the second term of equation 4.8 must be zero 
since Bessel functions of the second kind cannot be present with 
^ finite at the origin. 

d. Since the selection of the coordinate system is purely arbitrary 
the 0 = 0 axis can be positioned so that the phase angle in 
equation 4.7 is zero with no loss of generality. 



^^cLachlan, Bessel Functions for Engineers - 2nd Ed., Oxford at 
the Clarendon Press (1955) 

^^McLachlan, op. clt. 
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Thus, the solution of equation 4.2 takes the form: 

. ■~r 



A(TX 



- Fm.Kz'.kY & ^ e?"n( COS (VY^ 6 ) Jlrn(KYT)£- ^-9 

The final condition to be satisfied is that the pressure must be zero 
at r » C\. Thus 

3^ ( k.Y ~ o 

Defining a quantity S VY\Y\ such that 

S ~ ^ 

then 

'JrriCS^ 4=. ~ ^ 

where |T\ is an Integer denoting the nth zero of a Bessel function of order 
rn ♦ Further, it can be shown^® that all roots of equation 4.11 are 
real; that is, Syti^v must be real. Thus 

vy _ S> Yr\T\ 



is real. Since 



then 



>0 






'2- 



4.12 



Ky- = K - K 



K 






Substituting the complex forms for ^nd gives: 

— Kj: 2. 2. K'fc ^ ^ 0i\ + K "t* Z A!^\ K 

Equating the real and imaginary parts of equation 4.13 gives: 

=■ i- K' - 

^ 2 , Ki = o^, K' 

Combining equations 4.14a and 4.14b and solving for K gives: 

K'= ±'\|k^" r(^)^- ' 

If is assumed negligable. 






4. 14a 
4.14b 

4.15 

4.16 



18 

A. Gary and G. B. Mathews, A Treatise on Bessel Functions, Mac> 
millan & Co. Ltd. St. Martin's Street, London, pp 79, 1931. 
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where the negative root is rejected since only positive frequency is being 
considered. As might be expected, equation 4.16 is the equation for the 



It is desirable to excite a mode of vibration in which the medium 
has no tangential components of velocity. For a medium having even a 
slight viscosity, such motion at the walls can introduce a much greater 
damping of free vibrations than fhat due to absorption of energy by the 
medium itself. Thus, only axially symmetric modes will be considered 
(i.e. m * 0) for the remainder of the analysis. For the lowest axially 
symmetric mode, SYt^r^= So| “ *(0.7655) and = 1. This corresponds to 
a pressure antinode at the midpoint of the axis of the resonator. 

If attenuation is present, its effect will be to lower the resonant 
frequency slightly as indicated by equation 4.15. 



wave number in a lossless medium. Substituting for and in equa- 



tion 4.16 gives: 




Therefore, "F(iyr\Yx » the resonant frequency of the cavity is 




4.17 
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5. Sound Absorption and Dispersion By the Medium. 

Investigation of the general mechanism of absorption by the medium 
will be based on a solution of the first*>order equations obtained in 
section 3. They were 

^ *^o (V-V) — O (3.1) 5.1 

-?o t[jfoCv(T'i)//3]] (v- v) - J(v^e 

— "PoCv^S— O (3.3) 5.2 

-Po (H) + g-Q ^ -fofef )'! vf-/J N V(vv) 

_ ^ (3.5) 5.3 

ir/J'7x (yKV) - O 

where the subscripts indicating first-order effects have been (xnitted. 

Differentiating equation 5.3 with respect to time and substituting 
from 5.1 for gives 

5.4 

-^Nv(v-|:^)+/^vx(vx^) = o 

Considering that, according to Helmholtz's theorem^^ any vector can be 
uniquely separated into the sum of a gradient of a scalar and the curl of 
a divergenceless vector, let 

5.5 



where 



X 

V-A 



- velocity potential 

- vector potential 

n 0 



Substituting equation 5.5 into 5.4 yields the following scalar and vector 
wave equations: 

-Poll A =• 0 5.6b 

It is interesting to note that equation 5.6b which governs first-order 
vorticity is not a true vector wave equation but a vector analogue of the 






-V 



C 






^ - J3 •?, 



^ e 






5.6a 



^^Morse and Feshback, Methods of Theoretical Physics, Vol. 1, McGraw- 
Hill Book Co., pp. 53, 1953. 
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typical homogeneous diffusion equation. Since no source terms are present 
in 5.6b, this implies that, if first-order vorticity exists, it must be 
generated at the boundaries. This lack of source terms is due to the 
technique utilized in generating the first-order equations. Based on the 
assumption that the equilibrium values of the state variables are con- 
stant and that the viscosity coefficients are functions of these state 
variables, then it is consistent to assume that V/Jo «= vXo 0 • There- 
fore, the source terms involving the gradients of the viscosity coeffi- 
cients appear as second-order rather than fitst-order sources. This section 
deals specifically with absorption by the medium; thus, solution of 
equation 3.6b and its implications will be reserved for section 6. 

Neglecting the effects of radiation and assuming a time dependence 
of the form £ , equations 5.2 and 5.6a become: 

i-awa>n] h- 4 -A CO © = O 5.7 



+ -^wfoCv6 = O 



5.8 



Solving 3.7 for 0 and substituting into 5.8, and introducing the short- 
hand notation: 



and 



M/J (jJ 

G » ^ frequency dependent term 

^ o ^0 








“ ' /J WCp ' 


where 


-fo 


gives 






1 4-A&(l + YH) 



, which is non-frequency sensitive 






Equation 5.9 is biquadratic in the complex propagation constant and it is 
convenient to rewrite it as 






5.10 



The standard quadratic formula can be used to determine the reciprocal 
squares of the propagation constants as 
-Cq^ 










7' 






5.11 
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For the frequency range under consideration with water as the medium, G 
and H are small numbers and equation 5, 11 can be expanded by power aeries 
methods to yield: 






K,= 



60 



I -4- A 



r 



I + 



Neglecting second^order terms gives: 






A 6J 




k, = 



K,-= (i + a) 



i- 



(T-0 



K ~1 (a) 
■foC? J Co”^ 



-I k2- 

-Pq,Cp col ^ 

. 2. 4 J 



A/WCpji.'foCo^ 



5.12a 



5.12b 



5.13a 



5.13b 



Referring again to equation 5.10, we see that the solution will take the 
form of two waves referred to as types 1 and 11 whose wave numbers are 
respectively Ki snd ♦ Examination of equations 5.13 indicates that 

the two solutions might also be classified as "cempressional" and "thermal" 
since the type 1 waves represent the solution to the ordinary wave equa* 
tion with a slight damping coefficient and the type 11 waves are generally 
the result of thermal conduction. 

From equation 5.13a the absorption and dispersion measures for the 
type 1 waves are 



O i Cot) + ■|-(A’‘^N/-f„CoY] 



The attenuation thus determined is the Kirchhoff approximation for low 
frequencies. ^ 

The wave length of the type II wave ( is very much 

shorter than that of the type I wave. Further, the type II waves are 
very rapidly attenuated at the frequencies under consideration with a 



decrement of 
is given by: 



j. X 



"SCTT 



and the velocity of propagation through the medium 



>6 Cp" 
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6. Wall Losses. 

To study the effects of wall losses on the system, consider the force 
equation derived in section 2 (2.9). The viscosity coefficients /J and \ 
are assumed constant, and effects of the extraneous force, F, are assumed 
negligible. Further, since only small amplitude vibrations are to be 
considered, all terms involving squares of the velocity components are 
neglected. Equation 2.8 then becomes 



^ 4.(\ t v) + 



Since it was shovm in section 3 that the thermal waves are rapidly atten* 
uated, the effect of heat conduction at the boundary will be neglected. 
Thus, before proceeding further in the analysis, a direct relationship be* 
tween pressure and velocity must be established . Flow of mass out of 
any volume element reduces the pressure in a compressible fluid. In an 
elastic fluid, if compression is small, the rate of change of pressure is 
proportional to the divergence of the velocity; that is 



= -lb V* V 



6.2 



B - bulk modulus of elasticity 

Taking the time derlvitive of 6.1 and introducing the expression given in 

gives : 



equation 6.2 for 
x'i.rT 



Again defining V Qs 



6.3 



v-v4^ + vxA 



6.4 



and substituting equation 6.4 into 6.3 gives the two equations: 

^ -/jvxv 

Assuming a time dependence of the form where <T— , 

6.5a and 6.5b take the form: 



6.5a 

6.5b 

equations 



^^orse and Feshbach, Methods of Theoretical Physics, Vol. 1, McGraw- 
Hill Book Co., pp. 162, 1953. 
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VxV^A K^^ = o 
(v^ +- K^) = O 



6.6a 
6 • 6b 



where 



h 

k' 



2. -- A-P 

^ 2 
- -fo: 



■0)+- acf(zaj t X) 

satisfies a scalar Uelmholz equation, the form of solution will 



Since <|) 

be the sane as equation 4.9, that is 

4>=S\nK‘/2i JmCKyV) 



crV 



6.7 



where 

K-i - Si- 

3^ V(.y. C\^ =• O 

Equation 6.6a takes the general form of a vector "wave" equation. Accord* 
21 

ing to Stratton the three independent vector solutions to 6.6a are given 
by: 

U='^V 6.8 

^ = VX a ^ 6.9 

N ^ V X M 6.10 

where 0\ is any constant vector of unit length and 'V is a scalar func- 
tion satisfying the equation: 



From 6.8, 6.9 and 6.10, 

vIa = V'TJ = o 
v» ”u - ^ 

Since, by definition, 'V'* ji\=: O • 6.9 and 6.10 must be the complete solution 
of 6.6a; that is, ' 

^ "Fa 4-TT 

21j. A. Stratton, Electromagnetic Theory, McGraw-Hill Book Co., 
pp 393, 1941. 
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fe-; 

WM: 




1 






Equation 6.11 is again a scalar Helmholz equation and thus Is 

given by: 

Considering only the zero mode (m s 0) 

o;ur) 

Substituting equation 6*12 into 6*9 and 6.10 and selecting ^ gives 

T\- VK'o.-a '4^ 



N 



_ j__ 
" h 



(v X h) 



-I- C— wu s\n 
n [ S'^oo^ 



To(av') ^ *=■2^ 

Taking the curl o£ 6.13 and 6.14 gives 

vx A- a jbg'_^'i b=s 

+ Jo(a-) 

From the synmetry of the problem: 

= Vr (v-j l_- 

Thus, the rotational component of the vector velocity has the form 

'7XA= ‘a-r^b;8(bg “bin b£ Ji(gv)3 

■<s>D>'-,g(.bgcbs j, (av)] 



6.14 



6.15 



6.16a 

6.16b 



t* 0^1 



6.17 



where 






4 

From an examination of equations 6.7 and 6.17, it can be seen that since 
K'j.”' b> , then ^ ^ l^y . Thus 3"6 "9^ O • 

The irrotational component of velocity is determined by taking the 
gradient of 6.7. Thus 

V (t =r o^y [E. v(^ ,,V (-* ^Y- -^''1 ^ 2: 3 ; ( \< xr ^ 

+ 'a-2,[EviXr \C:.^ CUvV^J 



6.18 
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A rational accounting o£ the actual wall losses is well beyond the 
scope of this discussion. Derivation of even a reasonable approximation 
for the pressure field existing at the boundary would be difficult and 
tedious. Even If these equations were known, the mechanism of acoustical 
energy transfer to, and the absorption by the compliant structure is not 
sufficiently understood to make application feasible. Thus, it Is nec> 
essary to assume a simple model, considering that the solutions obtained 
will act as a guide In interpreting the actual behavior of the resonator. 

As a first approximation, the pressure release boundary condition 
Is retained with the additional assumption that the cylinder Is bounded 
by a smooth wall. It Is further assumed that, for a viscous medium, there 
Is no slipping between the medium and the wall. This implies the condl* 
tlon that the tangential components of velocity will be zero at the 
boundary . 

Examination of equation 6.18 Indicates that the Irrotatlonal com* 
ponents already satisfy this condition, l.e., 

Vr’O oX 

V 2 = 0 r = CK • 

Considering the rotational components, since 0, must be zero 

so that 

0 

Further, ^b,«i must equal zero for the longitudinal component, V 3 i to 
satisfy the boundary condition. Thus, for this simple model, there is no 
rotational velocity component, and no first-order vorticity exists. There- 
fore, the solution obtained In section 4 is unaltered. 

A closer approximation to the actual conditions can be obtained by 
assuming that pressure approaches some small finite value rather than zero 
at the boundary (l.e. , J“o( 0 ~0 ). In this case, the 

coefficient X)l 0 )g in equation 6.17 must again be equal to zero. However, 
from the two conditions: 



^ V-JL 


= 0 




V' ■= 


^ Vr 


= 0 


at. 


0 



two equations are obtained: 
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E. Kt S\' n ^ 3i^(KyV-)^ + O (8'^^= O 



6.19 



at 2-0 and 2 = L- ^ and 

E [Kz COS J' oC'^v^')'] + C^S^cosbi 3o( §c\) ~o. e. 

Since equation 6.19 must hold for all values of V • follows that 
Equation 6.20 can then be rewritten as 

To(v<t^v)- % T ^ 

~ ^ 3"6 • 
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kz-b-^^ 



6.21 



Equation 6.21 again indicates that for perfect pressure release walls (i.e. 

-0 ) there is no first»order vortlcity and no viscous loss of 
energy at the boundary. 

I 

As the boundary deviates from the Ideal case, the coefficient O in 
equation 6.21 Increases accordingly and vortlcity is generated. Thus, 
friction arises between the layer of fluid adhering to the walls and the 
iosnediately adjoining liquid in motion causing a dissipation of energy at 
the walls. 

This effect can be evaluated by considering that 



f 



Ky- ^ VY\V\ E VV\V^ 

where Ew\y^ is a small complex quantity which is dependent on O . Re- 
ferring to equation 4.13, it can be seen that the addition of £. will 
result in a decrease in the resonant frequency of the cavity as well as 
an increase in the total energy dissipation. 

A further Implication of this analysis is that, since O at T— (X * 
a transmitted as well as a reflected wave exists at the boundary and some 
energy will be lost to the surrounding medium due to radiation. 
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7. Measurement of Attenuation. 

Before proceeding to apply the results obtained in the preceding sec- 
tions to the actual design of a cavity, some method of measuring the 
energy absorption by the enclosed medium roust be determined. The quality 
factor or Q of a resonant cavity may be defined as: 



Q _ 2n (energy stored in cavity) 
“ energy lost per cycle 



7.1 



The mechanical energy stored in an incremental volume can be considered 
as the sum of the kinetic and potential energies; that is, 



- E 



K 



where 



- ±2 



^ Ep 
d V 



7.2 



^ ^ 'iL. 

Therefore^ W, the total energy stored in the cavity, is determined by tak- 
ing the volume integral of equation 7.2; i.e., 

-2. 



w= r c) V "E r c\ V 



7.3 



For the present, wall losses are neglected and all attenuation is assumed 

_ A<n> 

to occur in the medium. Selecting a time dependence of the form t, where 
CT is again equal to and substituting into equation 7.3 Indi- 

cates that W decays at a rate determined .by 0, where is the atten- 

uation constant in nepers per second. Thus, the energy lost during one 
period is given by: 

3.,CX W 



energy lost per cycle = 






7.4 



Substituting equations 7.3 and 7.4 into 7.1 gives 
n 2 -Trf W _ 

^ “ 2. W 2_ C>C 



or 












nepers/second 



7.5 



To convert 0^ from a time to a space dependence, equation 7.5 is multi- 
plied by : 

' k' 






lcQ ~ '2.G^ 



nepers/unit length 



7.6 



where K is defined by equation 4.15 as 
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Substituting In equation 7.6 and solving for OC gives: 






K: 



-t 



(^y 



7.7 



Thus, with the assumption that all energy dissipation occurs in the medium, 
the quality factor of the cavity can be utilized as a measure of the at- 
tenuation constant of the medium. 

In the case of an actual resonator, the effect of wall loss cannot 
be neglected, and the cavity is characterized by a total Q, Q^, given by 



J_. J_ * _L_ 

Q O. vy 

where Qw^ is due to dissipation by the medium and Q vv is a correction 
due to wall loss. Qm * then, is the value which must be utilized in 
equation 7.7 to measure absorption by the medium. 

Since the device under consideration will be operating at relatively 
low frequencies, determination of the correction factor, > is greatly 
simplified. 

It is well known that, at frequencies below 10 k.c./sec, attenuation 
of sound waves in fresh water is considerably less than attenuation of 
sound in sea water. In fact, by comparison, attenuation by the fresh 
water can be neglected. Thus, if the cavity is immersed in distilled 
water and the total quality factor is determined, then 

Qtcj Qwa 

to a close approximation. The subscript is used to indicate measure- 
ment in distilled water. 

Q may also be defined as 




where 



- resonant frequency 

- -T, • 

* frequencies of the upper and lower half power points. 



Therefore, 



Q 






-f )t 



A 






-f'cJ 



w 



Considering that the resonant frequencies in fresh and salt water will be 



only slightly different. 
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where 



- A-Cd 




Af 



band width between half power points for actual cavity 
in sea water 




band width between ha If -power points in distilled water 



f 



resonaitt frequency in sea water. 



Substituting this value into equation 7.7 gives 



TT ( \ 






yi-j H-n-cx) 



7.8 



and equation 7.8 can be used to experimentally determine attenuation by 
the medium. 
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8. Cavity Design. 

In the design and construction o£ a sound 'velocity and attenuation 
meter based on the principles discussed In the foregoing analysis, two 
general problems nnist be considered: 

a. Design of a cylindrical cavity presenting essentially a pressure- 
release boundary over the frequency range under consideration. 

b. Development of the "Intelligence" system Including a transducer 
to excite the water column to resonance, a receiving transducer 
to sample the generated pressure field and the associated 
electronic equipment for Interpretation and display of the de- 
sired Information. 

As a first step In the design, the general usage of the data obtained 
from the Instrument must be considered to determine the desired resonant 
frequency of the cavity. Equation 4.17 can then be utilized to calculate 
the approximate dimensions of the cylinder. Care must be taken to select 
a length to diameter ratio that will minimize danger of mode skipping. 

Next, the general range of sound velocities to be measured must be 
specified. This, In turn, determines how broad the Q of the cavity must 
be to present adequate sensitivity at the extremes of the operating range. 

With these parameters In hand, the formulas developed by Toulls^^ 
can be utilized to determine the nature and spacing of the compliant 
elements for construction of the cavity walls. This, then, outlines the 
steps to be carried out In the solution of the first portion of the 
design problem. 

In selecting and mounting the transducers for exciting the water 
column Into resonance and sampling the generated pressure field, several 
factors must be considered. First of all, since the material generally 
used for enclosing the transducers Is highly absorptive. It Is desirable 
to mount them In a relatively low Intensity sound field. Further, It 
Is necessary to keep separation sufficient to minimize direct coupling. 

The design of a system for Interpreting and presenting the data ob- 
tained Is primarily dependent upon the actual operating requirements and 
will not be discussed. 

J. Toulls, Acoustic Refraction and Scattering With Compliant 
Elements, to be published. 
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9. Conclusions. 

Exjualnation of the foregoing analysis Indicates that the solution 
to the simple wave equation obtained in Section 4 describes, to a very 
close approximation, the conditions existing in an actual resonator. 

From Section 5, it can be seen that the effects of heat conduction on 
the resonant frequency and attenuation constant are negligible. Section 
7 indicates that, with walls closely approximating an ideal pressure- 
release surface, the effects of first-order vortlclty can also be neg- 
lected. Thus, for design of a system employing the techniques under 
discussion, the simplified analysis presented in Section 4 should prove 
adequate . 

One particularly desirable feature of this system is that it mea- 
sures velocity and attenuation simultaneously. This should be of great 
assistance In determining the effect of air bubbles and other natural 
phenomena on sound velocity. 

Another desirable feature of this equipment is the frequency range 
employed. Since some doubt exists as to the validity of adapting sound 
velocity data measured at megacycle frequencies to much lower-frequency 
applications, an accurate device operating at these lower frequencies 
should be a welcome addition. 
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APPENDIX I 



CAVITY CONSTRUCTION AND TEST 

A cylindrical acoustic resonator sound>veloclty laeter was constructed 
by Toulls and preliminary tests were conducted by Toulls and the author 
from 6 January 1957 to 15 March 1957. This appendix outlines the methods 
of construction and Instrumentation utilized and describes the test 
results . 

The pres8ure*release walls of the cavity under test were constructed 
by deforming commercially available aluminum tubing (Alcoa 3003*H»14, 1/2 
hard, 1" O.D., 0.036" wall thickness) Into nearly flat or elliptic cross* 
section and welding the ends to prevent flooding. These tubes were then 
secured to a cylindrical frame with 2*1/2 Inches spacing between tubes 
forming a closed cylinder. 

The water column was excited Into resonance by a cylindrical barium 
tltanate transducer mounted Just below the compliant tube top of the 
cavity near the cylinder wall, and a similar transducer utilized to detect 
the pressure field was mounted directly opposite. 

For making velocity measurements, the output of the receiving trans* 
ducer was amplified and applied to the exciting transducer. A self* 
excited oscillating loop was thus formed. The loop voltage was applied, 
through a frequency multiplier (XIO), to a Potter 150 KC. Frequency 
Counter . 

For attenuation measurements, the amplitude of oscillation Is an 
Indirect measure of Q. Direct determination was made by driving the cavity 
with a General Radio Beat Frequency Oscillator (Type 1304*A) and roeasur* 

Ing the bandwidth between the half*power points. 

During the tests, the output of the receiving transducer was applied 
to the amplifier through a calibrated potentiometer. The mpllfler out* 
put was monitored and the gain (potentiometer) setting required to main* 
tain constant output was recorded as an Indication of the attenuation. A 
block diagram of the system Is shown In Figure 1. 

To evaluate the ability of the Instrument to measure actual sound 
velocity at sea, a number 'of tests were conducted In the coastal waters 
of San Diego, California. For comparison jpurposes and, particularly, to 
Isolate the effects of aallnlty, an additional set of runs was made at 
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the NEL testing facility at Sweetwater Lake. At each coastal station, 
the cavity was lo^jered by means of a bathythermograph winch, and the 
following data was recorded as the chamber was stopped at Intervals on 
the ascent: 

a. Resonant frequency as indicated by the Potter counter. 

b. Temperature as determined by a thermistor bead attached to the 
side of the cavity, 

c. Gain setting. 

d. Depth as indicated by a wire metering device. 

The results of these tests^^ Indicate that the system tested pre- 
sents promise for development as a survey Instrument for making sound- 
velocity profiles. 
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J. E, Wllle, Preliminary Test of A Cylindrical Acoustic Resonator 
Sound-Velocity lieter, NEL Tech. Memo., to be published. 
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